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CALCULATIONS  ON  DISSOCIATIVE  ATTACHMENT  TO  H-CL. 


Abstract • 

*  Dissociative  attachment  H  Cl  +  e  - >  H  +  Cl-  is  shown  to  occur  with 

cross-sections  of  order  1  angstrom~2  or  larger  only  at  impact  energies  within 
less  than  about  1  ev  of  the  threshold.  Cross-sections  for  DA  are  calculated 
fo-r  initial  vibrational  and  rotational  states  populated  significantly  at 
temperatures  up  to  1200  deg  K.  The  calculated  cross-sections  reproduce  the 
observed  result  that  raising  the  initial  vibrational  level  from  v-0  to  v-2 
raises  the  cross-section  by  a  factor  of  several  hundred*  The  energy  dependence 
of  the  c r os s-se c t i on  for  DA  to  molecules  at  temperatures  up  to  1200  deg  K 
agrees  in  orders  of  magnitude  with  experiment  for  both  HC1  and  DC1. 


CALCULATIONS  ON  DISSOCIATIVE  ATTACHMENT  OF  ELECTRONS  TO  H-Cl 


#1.  Introduction. 

The  object  of  the  contract  was  to  calculate  the  cross-section  for  the 

reaction  e  +  H  Cl  - >  H  +  Cl-  at  impact  energies  below  5eV.  The  effect  of  an 

initial  vibrational  excitation  of  the  molecule  was  to  be  studied;  initial 
levels  between  v®0  and  v**5  have  been  considered.  As  a  byproduct,  we  have 
calculated  vibrational  excitation  cross-sections  for  molecules  in  different 
initial  states  . 

The  physical  principles  and  the  method  of  calculation  are  described  in 
#2.  The  theory  is  tested  by  comparison  with  the  available  experimental  data 
in  #3.  The  remainder  of  the  calculated  cross-sections  is  given  in  #4.  #5  is 

a  short  summary. 

Some  important  magnitudes: 

Vibrational  quantum  in  HC1«0.371  ev ,  in  DC1*0.259  ev . 

Dissociation  threshold  of  cpld  HC1«0.81  ev ,  of  cold  DC1*0.87  ev . 

#2.  Physical  Principles  and  method  of  calculation. 

#2.1  The  model. 

The  observed  cross-sections  for  dissociative  attachment  to  the  molecule 
in  the  state  v*0  are  shown  in  Fig.l  (ref.l).  Evidently  there  is  only  a  single 
peak  at  impact  energies  below  5  eV .  Fig. 2  shows  the  results  of  published  ab 
initio  calculations  of  the  states  of  H  Cl  and  H  Cl-  (refs.2,3);  it  appears 
that  the  peak  in  Fig.l  between  1  and  2  eV  has  to  be  associated  with  a  single 
state  which  dissociates  into  H  +  C1-.  We  have  assumed  that  this  state  is  also 
the  one  responsible  for  the  peak  in  the  observed  cross-section  for  vibrational 
excitation  below  1,  which  is  shown  in  Fig. 3a  (ref. 4). 

The  calculations  done  under  the  contract  were  based  on  a  model  which  had 
previously  led  to  the  results  marked  'calculated'  in  Fig. 3b  (ref. 5).  The 
model  is  defined  by  the  following  points  (i)-(vi): 

(i)The  electron  is  treated  as  free  outside  the  molecule.  The  interaction 
with  the  tail  of  the  dipole  potential  is  neglected  because  this  interaction 
had  only  a  negligible  effect  in  the  calculation  of  vibrational  excitation  in 
r ef . 5  . 


(ii)Inside  the  molecule,  the  electron  is  fast  compared  with  the  nuclei, 
even  when  the  kinetic  energy  of  the  electron  outside  is  only  a  few  millivolts; 
therefore  the  wavefunction  of  the  electron  follows  the  nuclei  adiabat ical ly 
inside  the  molecule,  so  that  the  logarithmic  derivative  at  the  molecular 
surface  is  a  function  of  the  positions  of  the  nuclei,  and  independent  of  the 
nuclear  velocities. 


(iii) Outside  the  molecule,  the  wavefunction  of  the  extra  electron  may  be 
expressed  as  a  superposition  of  eigenstates  of  the  angular  momentum;  the 
different  1-values  are  coupled  at  the  molecular  surface  by  the  asymmetry  of 
the  potential.  At  energies  below  5  eV,  the  radial  vavef unctions  associated 
with  the  different  l-values  are  small  close  to  the  molecule  because  of  the 
centrifugal  potential  barriers.  The  value  1*0  is  an  exception,  so  that  its 
contribution  to  the  scattering  is  dominant  at  energies  below  5  eV.  Only  the 
s-wave  is  taken  into  account  when  the  electron  is  outside  the  molecule.  This 
approximation  is  consistent  with  the  spherical  symmetry  of  the  observed 
angular  distribution  of  scattered  electrons  in  vibrational  excitation  below  1 
eV  (ref .4) . 

(iv) The  energies  of  interest  are  sufficiently  low  for  the  logarithmic 

derivative  of  the  radial  function  in  the  s-wave  at  the  molecular  surface  to  be 

treated  as  independent  of  energy,  because  the  varying  kinetic  energy  of  the 
electron  outside  the  molecule  is  much  smaller  than  the  depth  of  the  potential 
inside.  The  logarithmic  derivative  of  the  s-wave  at  the  molecular  surface  may 
be  written  fO+fl(R-RO),  where  R  is  the  nuclear  separation,  and  RO  the  value  of 
R  at  equilibrium.  fO  and  fl  are  constants,  determined  from  a  fit  of  the 
observed  cross-sections  for  vibrational  excitation  (ref. 5).  The  linear 
approximation  with  respect  to  (R-RO)  is  justified  for  vibrational  excitation, 
where  the  excursion  of  the  nuclei  from  equilibrium  is  much  smaller  than  one 
bohr  (*0.52  angstrom)  during  the  collision,  because  in  the  absence  of  a 
trapping  mechanism,  the  collision  does  not  involve  any  time  delay.  (In  point 

(v)  below,  we  shall  return  to  the  range  of  validity  of  the  linear 

approximation.)  fl  is  of  order  unity  because  the  interior  electronic 
wavefunction  changes  by  a  substantial  fraction  of  itself  when  R,the  nuclear 
separation,  changes  by  one  bohr;  we  have  used  the  value  0.45  in  our 
calculations.  fO  is  close  to  zero,  with  a  fitted  value  of  0.075;  (this  value 
is  appropriate  if  the  logarithmic  derivative  refers  to  (r  x  6-wave),  rather 
than  to  the  s-wave  without  the  factor  r,  as  in  ref. 5). 

The  small  value  of  fO  implies  that  the  s-wave  of  the  electron  for  nuclei 
fixed  at  their  equilibrium  positions  is  approximately  r [ c os ( kr ) / ( kr ) J ,  instead 
of  r [ s in ( kr ) / ( kr  )  J  as  in  the  incident  plane  wave.  Thus  the  wavefunction  of 
the  electron  close  to  the  molecule  is  greatly  enhanced  relative  to  the 
incident  plane  wave;  this  enhancement  is  responsible  for  the  very  large 
cross-sect  ion  for  vibrational  excitation  near  threshold  (see  Fig. 3). 

(v) The  potential  for  the  negative  ion  can  be  calculated  from  the 
parameters  fO  and  fl,  the  binding  energy  being  ( f 0+f 1 ( R-RO ) ) "2/ 2 ;  the  result 
is  shown  in  Fig. 4.  The  potential  corresponding  to  the  fitted  logarithmic 
derivative  is  the  curve  marked  H  C1-,  which  leaves  the  curve  H  Cl  at  Rs*2.58 
bohr,  and  runs  to  the  point  A.  The  fit  of  the  vibrational  excitation  is 
sensitive  to  the  potential  only  close  to  the  equilibrium  separation;  however, 
the  linear  approximation  must  be  quite  good  out  to  A  because  the  fit  of  the 
transition  vm0  -->  v*2  is  satisfactory  in  ref. 5  ;  the  v*2  vibrational  level 
has  its  turning-point  0.45  bohr  to  the  right  of  the  equilibrium. 


There  cannot  be  a  barrier  beyond  A,  because  dissociative  attachment  is 
observed  to  set  in  strongly  at  threshold  (ref.l);  therefore  the  potential  must 
turn  over  near  A  to  a  nearly  horizontal  curve,  as  shown  in  Fig.4.  The 
resultant  potential,  represented  by  the  complete  curve  marked  H  Cl-  in  Fig.4, 
is  not  very  different  from  the  calculated  H  Cl-  potentials  given  in  refs. 2  and 
3. 


(vi)  For  dissociative  attachment,  ve  assumed  that  the  transition  to  the 
bound  H  Cl-  state  is  driven  by  the  non-adiabat ic  coupling  of  the  extra 
electron  to  the  velocities  of  the  nuclei  (Fig. 5).  In  other  words,  the 
electron  is  captured  into  the  bound  state  to  the  right  of  Rs  in  Fig. 5,  by 
transferring  its  binding  energy  to  the  nuclei.  This  mechanism  is  quite 
different  from  the  familiar  mechanism  of  dissociative  attachment  at  a 
quasistat lonary  state,  where  the  potential  curve  for  the  negative  ion 
continues  from  the  region  of  finite  binding  for  the  extra  electron  into  the 
region  of  positive  energies;  here  attachment  starts  when  the  electron  tunnells 
into  the  molecule  in  a  manner  which  does  not  depend  on  the  nuclear  velocities. 
In  the  absence  of  a  trapping  mechanism  at  positive  energies  for  the  extra 
electron,  there  is  no  localised  electronic  state  which  one  might  associate  an 
H  Cl-  ion  when  R<Rs. 

#2.2  Method  of  calculation. 


A  wavefunction  which  includes  dissociative  attachment  may  be  written 

Tp'iVR’)  -  %C£)  +  U rC<J'£)f  (1 

where  Cf  is  the  bound  H  Cl-  eigenfunction  of  the  electronic 
Hamiltonian;  it  belongs  to  the  eigenvalue  is  the  associated 

nuclear  wavefunction  of  the  negative  ion.  One  may  write 


VR>  -  Ccf  l2) 

where  the  notat  ion£...  U..)^means  that  the  integral  extends  over  the 
electron  co-ordinates  while  the  nuclei  are  held  fixed,  is 

the  part  of  the  wavefunction  orthogonal  to  Cf  in  the  space  of  electron 
co-ordinates.  One  has  because  ®  at  finite  distancesTof  the 

electron  when  R*Rs  in  Fig.4,  due  to  the  infinite  radial  extent  when  the 
binding  energy  vanishes. 


If  one  inserts  into  the  Schroedinger  equation,  multiplies  by  CJf . , 

ana  integrates  over  the  co-ordinates  of  the  electrons,  one  obtains  the  nuclear 
wave-equation  for  the  negative  ion: 


\  X5 

2KI 


There  cannot  be  a  barrier  beyond  A,  because  dissociative  attachment  is 
observed  to  set  in  strongly  at  threshold  (ref.l);  therefore  the  potential  must 
turn  over  near  A  to  a  nearly  horirontal  curve,  as  shown  in  Fig.4.  The 
resultant  potential,  represented  by  the  complete  curve  marked  H  Cl-  in  Fig. 4, 
is  not  very  different  from  the  calculated  H  Cl-  potentials  given  in  refs. 2  and 
3. 


(vi)  For  dissociative  attachment,  we  assumed  that  the  transition  to  the 
bound  H  Cl-  state  is  driven  by  the  non-adiabat ic  coupling  of  the  extra 
electron  to  the  velocities  of  the  nuclei  (Fig. 5).  In  other  words,  the 
electron  is  captured  into  the  bound  state  to  the  right  of  Rs  in  Fig. 5,  by 
transferring  its  binding  energy  to  the  nuclei.  This  mechanism  is  quite 
different  from  the  familiar  mechanism  of  dissociative  attachment  at  a 
qua s is t at xonary  state,  where  the  potential  curve  for  the  negative  ion 
continues  from  the  region  of  finite  binding  for  the  extra  electron  into  the 
region  of  positive  energies;  here  attachment  starts  when  the  electron  tunnells 
into  the  molecule  in  a  manner  which  does  not  depend  on  the  nuclear  velocities. 
In  the  absence  of  a  trapping  mechanism  at  positive  energies  for  the  extra 
electron,  there  is  no  localised  electronic  state  which  one  might  associate  an 
H  Cl-  ion  when  R<Rs. 

#2.2  Method  of  calculation. 


A  wavefunction  which  includes  dissociative  attachment  may  be  written 

Tp-<VR')  -  TC£)  +  t-'c-o-r),  (1> 

where  Cp  is  the  bound  H  Cl-  eigenfunction  of  the  electronic 
Hamiltonian;  it  belongs  to  the  eigenvalue  is  the  associated 

nuclear  wavefunction  of  the  negative  ion.  One  may  write 

VR)  *  Ccf  <2) 

where  the  notation^.-  I.M^means  that  the  integral  extjends  over  the 
electron  co-ordinates  while  the  nuclei  are  held  fixed.  is 

the  part  of  the  wavefunction  orthogonal  to  Cf  in  the  space  of  electron 
co-ordinates.  One  has  because  0  at  all  finite  dist ancesTbf  the 

electron  when  R»Rs  in  Fig. 4,  due  to  the  infinite  radial  extent  when  the 
binding  energy  vanishes. 


If  one  inserts  into  the  Schroedinger  equation,  multiplies  by  Cft  • , 

and  integrates  over  the  co-ordinates  of  the  electrons,  one  obtains  the  nuclear 
wave-equation  for  the  negative  ion: 


-4- 


i s  tbe  reduced  mass  of  the  nuclei.  One  has 


x  because  C\f\ *^0£  *1,  because  of  the  normalisation.  The  non-adiabat ic  term 


C'fl  is  small  compared  with  the  kinetic  energy  term 

because  c^>  varies  slowly  with  R  in  comparison  with  c  **** 
noh-adiabat ic  term  was  dropped.  5 


this 


The  amplitude  of  ^  when  R  - >  oO  is  given  by 

1  — > 


^  ^  ^  ^  1  ?>&  - 

0  K* 

The  function  is  the  solution  of  the  homogeneous  version  of  the  nuclear 

wave-equation,  which  satisfies  *0,  and  behaves  asymptotically  like 

- >  6<*^  )  .  K  is  given  by  K*  (  £“o6*i)  -  £~] ,  ^ «  The 

derivation  uses  the  boundary  conditions  ^  - as  R - >  cx^  , 

and  (r,R, )  *0.  .In  the  calculations  under  the  contract, 

the  integral  fctRT.  .  .  was  worked  out  with  the  asymptotic  form  for  the 
electronic  wavefunction  qp  ,  because  most  of  the  transition  must  occur  at 

small  binding  energies,  where  the  radial  extent  of  the  electronic  wavefunction 
is  very  large  compared  with  the  target.  The  second  order  term  containing 


5'cf/3  ^ 


was  dropped. 


For  the  function  ,  we  used  the  wavefunction  determined  by  the 

solution  of  the  boundary  condition 


at  the  molecular  surface  r£ 


by  an  expansion  of  the  wavefunction  in 


vibrational  eigenfunctions  of  the  neutral  molecule.  (That  expansion  yields 
the  vibrational  excitation.)  Nine  vibrational  levels  were  taken  into  account; 
this  number  was  chosen  after  some  tests  to  check  on  convergence.  The  boundary 
condition  was  applied  at  r ^m0  ,  because  the  wavefunction  of  the  electron 
outside  the  molecule  varies  on  a  scale  large  compared  with  the  size  of  the 
molecule. 

,  A  peculiar  difficulty  arises  in  the  calculation  of  this  approximation  for 
Because  the  tails  of  the  vibrational  wavef unctions  in  the  neutral  stick 
out  into  the  region  in  Fig. 4  where  the  HC1-  potential  lies  below  the  HC1 
potential,  there  are  quasistationary  states  due  to  the  attachment  of  an 
electron  to  the  excited  vibrational  levels;  6uch  states  are  sometimes  called 
'nuclear  excited  Feshbach  resonances'.  The  resonances  of  this  kind  which  are 
associated  with  the  levels  v*l  and  v*2  are  probably  genuine  physical  phenomena 
(ref. 6);  they  lie  below  the  threshold  for  dissociation  into  H  +  Cl-  ,  so  that 
the  emission  of  an  electron  is  their  only  mode  of  decay.  However,  if  one  uses 
the  potential  for  HC1-  associated  with  the  linear  approximation  for  the 
logarithmic  derivative,  then  the  part  of  the  HC1  potential  which  extends 
beyond  the  point  A  in  Fig. 4  gives  rise  to  quasistationary  states  due  to  the 
attachment  of  an  electron  to  the  HC1  levels  from  v*3  upwards,  in  addition  to 


the  levels  below  v*3  .  But  the  higher  levels  are  artifacts,  because  the 
approximation  for  the  HC1-  potential  from  the  linear  logarithmic  derivative 
does  not  allow  them  to  decay  by  dissociation^,  ^  which  is  now  energetically 
possible.  If  one  calculates  the  function  in  the  manner  of  ref.5,  with 

the  linear  approximation  for  the  logarithmic  derivative,  then  it  resonates  at 
the  artificial  resonances.  The  artificial  resonances  lie  too  high  to  affect 
thfe  vibrational  excitation  in  the  threshold  region  discussed  in  ref.5,  but 
they  introduce  very  serious  errors  in  the  dissociative  attachment  which  occurs 
above  the  dissociation  threshold. 


A  really  satisfactory  treatment  of  the  artificial  resonances  was  not 
possible  within  the  limited  scope  of  this  contract.  What  is  needed  is  a 
calculation  of  the  vibrational  excitation  and  dissociative  attachment 
simultaneously,  with  a  final-state  potential  for  HC1-  which  reproduces  the 
dissociation  energy.  In  such  a  treatment,  the  artificial  resonances  would  be 
destroyed  by  dissociation.  However,  that  would  be  a  much  more  ambitious 
calculation  than  the  treatment  of  dissociative  attachment  as  a  perturbation  on 
vibrational  excitation,  which  is  what  I  had  planned  for  this  contract.  I 
intend  to  continue  working  on  the  more  ambitious  theory,  but  for  the  present, 

I  have  used  a  'crutch'  to  get  rid  of  the  artificial  resonances.  The  crutch 
retains  the  linear  approximation  for  the  logarithmic  derivative  to  calculate 
tne  function  Tfr*  ,  but  drops  the  contribution  of  any  vibrational  state  of  HC1 
from  the  function  whenever  the  total  energy  is  not  more  than  100  mev 

below  that  state.  This  trick  never  drops  more  than  a  few  of  the  levels 
contributing  to  the,  function  "Ip**  ,  because  the  spacing  of  the  vibrational 

levels  is  350  meV,  and  because  the  different  rotational  sublevels  of  any 
vibrational  level  differ  by  up  to  about  200  meV  in  rotational  energy.  The 
accuracy  of  the  results  can  be  judged  by  comparison  with  the  limited  available 
measurements;  (see  #4). 


#3.  Tests  of  the  theory:  Comparison  with  observations* 

#3.1  Vibrational  Excitation. 

*  Our  calculated  results  for  vibrational  excitation  near  threshold  are 
shown  in  Fig. 3b.  Fifteen  vibrational  levels  of  BC1  were  included  in  the 
calculation.  (The  calculations  on  DA  used  9.)  There  is  nothing  novel  about 
th.ese  results;  they  merely  show  that  the  programs  used  for  the  present 
calculations  were  at  least  as  good  at  calculating  vibrational  excitation  as 
t de  one's  Dube  and  I  used  four  years  ago  on  the  same  problem  in  ref. 5. 

The  calculation  uses  a  harmonic  oscillator  model. 

#3.2  Dissociative  Attachment. 

The  calculation  of  dissociative  attachment  sketched  in  the  previous 
paragraph  uses  the  wavefunction  behind  the  results  in  Fig. 3  to  evaluate  the 
integral  in  eq.5.  In  addition  to  this  wavefunction,  one  needs  a  potential  for 
the  nuclei  in  the  negative  ion.  We  used  the  potential  in  Fig. 4;  at 
separations  smaller  than  the  point  A,  it  follows  from  the  linear  approximation 
for  the  logarithmic  derivative  in  equ.6.  The  point  A  was  chosen  so  that  the 
potential  is  0.05  ev  below  the  asymptotic  limit.  From  A,  the  potential  was 
joined  smoothly  to  the  polarisation  potential  of  an  H  atom  near  the  single 
charge  of  the  Cl-  ion  at  large  separations.  The  rest  of  the  calculation 
involves  no  further1 adjustments  of  physical  parameters. 

The  calculation  was  done  for  different  initial  angular  momenta  and 
vibrational  states  of  the  nuclei,  up  to  1*18  and  v*2  for  HC1  and  up  to  1*27 
ana  v*4  for  DC1.  (As  we  said  before,  9  vibrational  levels  were  included  in 
the  calculation.  The  limiting  values  were  chosen  to  reach  all  the  angular 
momentum  states  populated  significantly  at  temperatures  up  to  1200  deg  K.) 

The  angular  momentum  of  the  nuclei  comes  entirely  from  the  initial  molecular 
state,  because  the  electron  is  presumed  to  enter  in  an  8-state. 

The  available  experimental  measurements  for  HC1  +  e  -->  H  +  Cl-  at 
temperatures  up  to  1200  degrees  K  are  shown  in  Fig. 6a,  taken  from  a  recent 
paper(7).  The  c or r e spona ing  calculated  results  are  shown  in  Fig. 6b;  these 
curves  were  obtained  by  a  Boltzmann  average  of  the  calculated  results  for 
different  initial  1  and  v  for  the  molecule.  The  electron  impact  energy  was 
varied  in  steps  of  50  mev;  each  cross-section  was  smeared  out  over  its  50  mev 
interval  in  the  pictures. 

It  is  important  that  any  initial  excitation  energy  of  the  molecule,  in 
vibration  or  rotation,  contributes  to  the  dissociation  energy.  Therefore 
the  dissociative  attachment  thresholds  for  the  excited  states  are  lower  than 
the  threshold  for  the  ground  state  by  the  initial  excitation  energy.  The 
threshold  in  HC1  is  0.81  ev  for  the  ground  state  (v*0,  1}0),  and  0.44  ev  and 
0.07  ev  for  (v*l,l*0)  and  (v*2,l*0).  The  difference  in  thresholds  makes  it 
easy  to  identify  the  contribution  of  excited  states  in  Fig.6  near  their 
thresholds.  Since  the  mean  rotational  kinetic  energy  is  only  0.1  ev  at  1000 
deg  K,  the  strong  cross-sections  in  Fig.6  well  below  the  cold  threshold  at 
0.81  ev  must  involve  initial  vibrational  excitation. 


1 


The  most  striking  aspect  of  the  results  in  Fig. 6  is  the  important 
contribution  of  the  excited  vibrational  states  of  the  target.  Because  of  the 
low  population  of  the  state  v*2  in  BC1,  the  dissociative  attachment 
cross-sections  to  the  state  v*2  must  be  a  few  hundred  times  larger  than  the 
cross-sections  for  v*0  to  show  as  large  in  the  experimental  curves  as  it  does. 
The  theory  reproduces  the  relative  importance  of  the  higher  v  levels  in 
reasonable  agreement  with  experiment.  There  is  also  an  indication  of  an 
enhancement  of  the  cross-sections  by  rotational  excitation  in  the  molecule, 
when  the  rotational  kinetic  energy  lowers  the  DA  threshold  through  one  of  the 
v lb r a t 1 ona 1 ly  excited  target  states.  This  enhancement  is  again  reproduced 
fairly  well  by  the  theory.  (The  peak  marked  '3'  on  the  calculated  curve  for 
1180  deg  K  is  due  to  this  effect.)  The  mechanism  of  this  effect  is  not 
completely  understood. 


Fig. 7  shows  the  observed  and  calculated  cross-sections  for  DA  in  DC1  at 
ll^u  deg  K  (observations  from  ref. 7).  There  is  again  a  very  large 
contribution  from  the  v ibra t iona 1 1  y  excited  states  of  the  target,  whose 
relative  magnitude  is  reproduced  by  the  calculation. 

The  experimental  measurements  in  Fig. 6a  and  Fig. 7  are  not  absolute.  The 
only  absolute  measurements  available  are  those  shown  in  Fig.l  ,  for  DA  to  cold 
HC1  and  DC1.  These  results  are  compared  with  the  calculations  in  Fig. 8,  where 
the  points  indicate  the  experiments,  and  the  continuous  curve  the  theory.  The 
calculation  lies  between  a  factor  2  and  3  above  the  experiments,  even  after 
one  allows  for  a  smoothing  out  of  the  structure  in  the  calculated  results. 

This  discrepancy  in1  the  absolute  magnitude  is  the  major  failure  of  the  theory. 
However,  the  ratio  of  the  cross-sections  for  HC1  and  DC1  in  Fig. 8  is  correct 
in  the  theory. 

#4.  Calculated  cross-sections  for  individual  channels. 

The  calculated  cross-sections  for  different  initial  and  final  channels 
are  shown  in  Fig. 9  for  vibrational  excitation,  and  in  Fig. 10  for  dissociative 
attachment.  These  are  the  cross-sections  which  are  combined  in  Figs. 3, 6, and  7 
with  Boltzmann  weighting  factors. 

Our  model  neglects  coupling  between  vibrations  and  rotations,  so  that  the 
rotations  enter  only  through  their  kinetic  energy. 

An  examination  of  the  calculated  cross-sections  shows  that  they 
do  not  satisfy  the  detailed  balance  theorem  exactly.  Generally,  and 
particularly  for  the  larger  cross-sections,  the  ratio  of  inverse 
cross-sect  ions  lies  within  a  factor  1.5  of  the  value  required  by  detailed 
balancing.  The  origin  of  this  difficulty  is  that  the  calculation  treats 
dissociative  attachment  as  a  perturbation  on  vibrational  excitation;  one  would 
have  to  treat  the  two  processes  in  a  more  symmetric  manner  to  satisfy  detailed 
balancing. 


#5*  Summary  and  conclusion. 

We  conclude: 

The  model  which  accounted  for  vibrational  excitation  of  HC1  in  the 
thresnold  region  in  ref. 5  accounts  also  for  the  following  aspects  of 
dissociative  attachment  HCl+e — >H  +  C1-: 

(1)  The  dependence  on  energy  ;  in  particular,  the  model  suggests  that  the 
production  of  Cl-  ions  should  occur  only  at  electron  impact  energies  within 
one  ev  of  the  threshold. 

(2)  The  ratios  of  the  cross-sections  for  different  initial  vibrational 
and  rotational  states;  in  particular,  the  model  predicts  an  increase  of  a 
factor  of  several  hundred  in  the  cross-section  when  the  molecule  is  raised 
from  the  vibrational  ground-state  to  the  second  vibrational  excited  state. 

(3)  The  ratios  of  the  cross-sections  for  H  Cl  +  e  -->  H  +  Cl-  and 
DC  1  +  e  — >  D  +  C1-;  in  particular,  the  model  suggests  that  for  the 

lowest  vibrational  state,  the  cross-section  should  be  about  five  times  smaller 
for  DC1  than  for  HC1. 

The  model  gives  DA  cross-sections  too  large  in  absolute  magnitude  by  a 
factor  of  2-3  in  the  only  case  where  measurements  have  been  made,  for 
molecules  with  no  vibrational  excitation. 

The  calculations  do  not  satisfy  detailed  balancing  exactly,  essentially 
because  they  treat  dissociative  attachment  as  a  perturbation  on  vibrational 
excitation.  For  the  larger  cross-sections,  the  ratio  of  the  cross-sections  for 
inverse  processes  generally  lies  within  a  factor  1.5  of  the  value  required 
by  detailed  balancing. 

The  calculations  show  that  the  DA  cross-sections  fall  rapidly  as  the 
energy  rises  above  threshold,  within  about  one  ev  or  a  little  less.  Together 
witn  the  fact  that  there  is  only  a  single  state  available  for  the  dissociation 
into  H  +  Cl—,  according  to  refs.2  and  3,  this  rapid  fall-  off  shows  that  the 
thresnold  peaks  of  dissociative  attachment  we  have  calculated  are  the  only 
contribution  to  dissociative  attachment  at  impact  energies  below  5  ev • 
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FIGURE  CAPTIONS 


Fig.l.  Observed  cross-sections  for  dissociative  attachment  to  H-Cl. 

(From  ref.l.)  The  cross-section  for  D-CI  +  e  — >  D  +  Cl-  is  similar  in 

shape  but  smaller  by  a  factor  5  at  the  threshold.' 

Fig. 2.  Ab  initio  calculated  energy- level s  for  H  Cl  and  H  C1-.  (From 
ref  8 .  2  and  3  .  ) 

Fig. 3a.  Dots  .‘Observed  cross-sections  for  vibrational  excitation 

e  +  H  Cl(v«0)  - >  e'  ♦  H  Cl(v*l).  Continuous  curve:  Calculated  in  1977. 

(From  ref. 5.  Data  from  ref. 4.) 

3b.  Calculation  with  the  present  program. 

Fig. 4.  Potential  energy  curves  used  in  the  present  calculation. 

H  Cl:  Harmonic  oscillator  model. 

H  C1-:  Binding  energy  of  the  extra  electron  from  the  logarithmic  derivative  of 
the  wavefunction  which  fits  the  vibrational  excitation  in  Fig. 3a,  at 
separations  smaller  than  A.  At  larger  distances,  polarisation  potential. 

Fig.5a.  The  traditional  adiabatic  model  for  dissociative  attachment, 
where  the  incoming  electron  tunnels  into  a  qua s i st a t ionary  state  through  a 
potential  barrier.  The  nuclei  then  slide  apart  with  the  extra  electron 
trapped  by  the  barrier,  until  the  electronic  state  becomes  truly  bound  at  the 
point  Rs. 

5b.  The  non-adiabat ic  model  behind  the  present  calculations.  There  is 
no  barrier  to  hold  the  electron  at  separations  smaller  than  Rs.  An  incoming 
electron  is  captured  into  the  bound  state  at  separations  larger  than  Rs  by  the 
coupling  to  the  finite  velocity  of  the  nuclei. 

Fig.6a.  Observed  DA  to  H  Cl  at  various  temperatures  up  to  1200  deg  K. 

(From  ref. 7.)  The  features  marked  in  the  top  diagram  may  correspond  to 
the  features  marked  in  Fig.6b. 

6b.  The  present  calculations  of  DA  at  the  temperatures  corresponding  to 
Fig. 6a. 

Fig. 7.  Observed  DA  to  D  Cl  (from  ref. 7),  and  corresponding  results  from 
the  present  calculation. 

Fig.8.  Absolute  magnitudes  in  experiment  and  theory.  Dots:  Observations, 
from  Fig.l  and  ref.l.  Continuous  curve:  Present  calculation. 

8a.  H  Cl. 

8b.  D  Cl . 

Fig. 9.  The  present  calculations  for  vibrational  excitation  from 
different  initial  vibrational  states,  all  for  1*0 . 

Fig.  10.  The  present  calculations  for  dissociative  attachment  for 
different  initial  vibrational  states. 

10a.  H  Cl. 

10b.  D  Cl. 
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